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We analyze the screening and bulk energy of a classical and strongly interacting plasma of color 
charges, a model we recently introduced for the description of a quark-gluon plasma at T — (1 — 3)T C . 
The partition function is organized around the Debye-Hiickel limit. The linear Debye-Hiickel limit 
is corrected by a virial expansion. For the pressure, the expansion is badly convergent even in 
the dilute limit. The non-linear Debye-Hiickel theory is studied numerically as an alternative for 
moderately strong plasmas. We use Debye theory of solid to extend the analysis to the crystal phase 
at very strong coupling. The analytical results for the bulk energy per particle compare well with 
the numerical results from molecular dynamics simulation for all couplings. 



I. INTRODUCTION 

QCD at high temperature is believed to be in a quark 
gluon plasma (QGP) phase, whereby color charges are 
screened rather than confined [1]. Asymptotic freedom 
guarantees that for T > A (the QCD cutoff) the QGP is 
weakly coupled (wQGP) with dressed quarks and gluons 
behaving as quasi-particles near the ideal gas limit. 

The QGP has been intensely sought by dedicated ex- 
periments at CERN SPS and more recently at BNL 
RHIC collider. Extensive analysis of RHIC data have 
revealed collective effects known as radial and elliptic 
flows. Their hydrodynamical explanation [3] have sug- 
gested that the QGP is a near-perfect liquid [4], promptly 
formed in heavy ion collisions at temperatures of the or- 
der of T w 2T C . As all dissipative lengths in the QGP 
seem to be short, it is not a weakly coupled gas-like phase 
but more like a liquid-like phase. 

Recently, two of us [5] have suggested that the in- 
teraction of quasiparticles in the relevant temperature 
range probed at RHIC is strong enough to form multiply 
marginal bound states. Some of these states (charmo- 
nium) were recently reported in lattice simulations [6] at 
temperatures as high as T sa 2T C . Many more colored 
and uncolored states were predicted, and remain to be 
analyzed by lattice simulations. The QGP in this regime 
will be referred to as a strongly coupled QGP (sQGP). 
First principle calculations of the sQGP properties can 
be done for supersymmetric extensions of QCD, N = 4 
SYM, via the AdS/CFT correspondence. 

In a recent paper [7], hereby referred to as I, we sug- 
gested a classical nonrelativistic model to help under- 
stand and describe some relevant features of the sQGP. 
Essentially it is the transport properties of the sQGP, 
which arc difficult or impossible to access through lattice 
simulations. The model consists of massive nonrelativis- 
tic quarks and gluons, which interact via color Coulomb 
interactions plus some repulsion, for stability. The color 
charges are assumed to be large and classical obeying 
Wong equations of motion. In a way, we are suggesting 
that quantum effects in the QGP are basically reduced 



to generate thermal-like masses, cause the effective cou- 
pling to run to larger values at smaller T, and add the 
"localization energy" to the Coulomb interaction. This 
model, referred to as classical QGP (cQGP) was studied 
in I using molecular dynamics simulations. 

In this paper we provide some analytical analysis of the 
static bulk properties such as the pressure and the energy 
density of the cQGP thereby unraveling its phase struc- 
ture at weak, intermediate and strong coupling, where is 
is in a gas, liquid and solid phases respectively. Similar 
approaches in the context of solutes of different natures 
where also considered [8]. 

In section II we review briefly the cQGP model. In 
section III we carry an exact low density expansion of 
the cQGP around the linearized Dcbyc-Huckcl limit to 
sixth-order in the particle concentration. In section IV 
we go beyond the linear Debye-Hiickel theory and study 
screening in a nonlinear regime. In section V we study 
another version of screening, via the crystal phase. In 
section VI we use Debye theory of solids to extend the 
Debye-Hiickel result for dilute but screened gases all the 
way to the crystalline phase. Our analytical results for 
the excess energy and pressure compare favorably to the 
results from molecular dynamics at all couplings. Our 
conclusions are summarized in section VII. 



II. CLASSICAL QGP MODEL 

For quarks and gluons with large thermal energies of 
the order of ttT, we may consider the cQGP as a system 
of non-relativistic charges interacting solely through lon- 
gitudinal color electric fields. The magnetic effects are 
subleading in the non-relativistic limit. Specifically, the 
Hamiltonian for the cQGP reads 
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where the sum is over the species a = q,q, g and their 
respective numbers N a , each carrying a thermal energy 
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m a i=a ttT. We have added a short-range repulsive core, 
which is needed to regulate the short-distance integrals 
and dynamical simulations *. The phase space coordi- 
nates are position (x a ), momentum (p a ) and color (Q a ). 
Only the Coulomb- like interaction was retained in (1) in 
the non-relativistic limit. The non-perturbative effects 
due to chromomagnetism will be discussed elsewhere. 

A central issue in classical colored plasmas is screen- 
ing. In the weak coupling limit the linear Debye-Huckel 
theory is well established and we will describe it below. 
At intermediate couplings (say T = 1 — 10 in the no- 
tations of I) not much is known and we will give some 
suggestions in this direction below. At large couplings 
(say r = 100 — 200 in the notations of I) a more ade- 
quate description can be developed starting instead from 
the crystal limit. 



III. CLASSICAL QGP PARTITION FUNCTION 

In this section we discuss the thermodynamics of 
cQGP. We formulate a generalized partition function for 
the cQGP to all orders in the classical coupling T (ra- 
tio of the potential to kinetic energy) and density. We 
discuss screening as a mean-field phenomenon. We first 
discuss its weak coupling limit in the form of a linear 
Dcbye-Hiickel limit, and corrections through a low den- 
sity virial expansion. 

The partition function of the cQGP is 
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and the density 
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The unscreened and dimensionless Coulomb potential is 
v(x) — Ib/t with Ib = g 2 /(AitT) the length at which two 
charges interact with energy T t. To prevent a classical 
collapse of the cQGP, we introduce a phenomenological 



* Minimally this term is caused by quantum localization 
energy. 

t We have substituted g 2 by p 2 /47r to make contact with 
standard SI conventions in strongly coupled plasmas 



repulsion with a core a that is w(x) = oo for x < a 
and zero otherwise. In a quantum theory a is set by the 
Heisenberg uncertainty principle. The A Q 's in (2) are just 
the thermal wavelengths for each species a = q,q,g. 

At weak coupling or in the gas phase, the partition 
function (2) is expected to scale as 
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The 3/2 is standard from the Maxwellian distribution 
over momenta, while a(N c = 2) = 1 and a(N c = 3) = 3 
are extra powers from the classical color degrees of free- 
dom. Indeed for iV c = 3 there are 3 angles and 3 con- 
jugate angles (generalized momenta) in addition to the 
2 fixed Casimirs that make the 8-dimensional classical 
color vector. 

At very strong coupling or in the crystal phase, the 
partition function (2) is expected to scale as 
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since the crystal localizes the particles in space as well. 
The change in behavior between (5) and (6) is reflected 
in the change in the specific heat of each phase. The 
partition function and thus the specific heat of the liquid 
is intermediate between (5) and (6). 



A. Mean-field screening 

The grand-partition function follows from (2) through 
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where the classical fugacities for the quarks, antiquarks 
and gluons are the same as the three species are all in the 
adjoint representation with comparable thermal energies 
m a . 

Since (2) is Gaussian in the densities it can be lin- 
earized using standard Hubbard-Stratonovich transfor- 
mations. The result after exponentiation is 
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with the induced action 



= \J dxdx'(f> a {x)v- 1 (x-x')(f> a (x') 
+ i J dxdx'ip(x) w^ 1 (x — x') tp(x') 
+A ^ dQ a dx a 

a 

x e i(w(0)+v(0)+Ql-iQ a a 4> a (x a )-ii(>(x a )) 



(9) 



The contributions v(0) and w(0) are the divergent self- 
energies. The normalizations Z ViW are 
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This analysis borrows on the functional framework de- 
veloped for various ionic systems as discussed in [8] (and 
references therein). 

The induced action (9) captures the essentials of the 
cQGP at all couplings. It is instructive to see that the 
saddle point approximation 5S/5<fi = with ip = yields 
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If we define <fi = (—ig/ksT) 0, then (11) reduces to 
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which is just the non-linear Debye-Hiickel equation for 
the cQGP with the renormalized concentration 
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For one species, a trivial mean-field solution is Cp = Qip 
with VV = 0. 

The mean-field equation (12) is valid for both weakly 
and strongly coupled plasma. At weak coupling it re- 
duces to the linear Debye-Hiickel theory which we will 
explore below. At intermediate and strong coupling it 
corresponds to the non-linear Debye-Hiickel description 
of which the crystalline structure is an asymptotic limit. 
This will also be discussed below. 



B. Linear Debye-Hiickel Screening 

The linearized Debye-Hiickel limit follows from (12) 
for weak coupling or high temperature, whereby the ex- 
ponent is expanded to first order leading to 
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This suggests to organize (9) around the Debye-Hiickel 
limit. Specifically, we introduce the screened Coulomb 
potential 
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with K 2 = 4irl B \J2 a Ql/(N 2 - 1). To leading order A 
is the concentration of tertiaries 
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as we will show below. So K 2 is just the sum of the second 
Casimir's weighted with the species concentration. Using 
the inverse relation 
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We may reorganize (8) around the linearized Debye- 
Hiickel limit. Specifically 
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The averaging in (18) is carried using the linearized 
Debye-Hiickel measure 
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We note that (18) while organized around the linearized 
Debye-Hiickel limit, it still includes the non- linear cor- 
rections in the coupling to all orders. 

In the linearized Debye-Hiickel limit, the dimcnsionlcss 
pair interaction is 

0(1, 2) w e - Q i Q 2 »dh(*i-Ro) (22) 

and fulfills the superposition principle in liquids 

.9(1,2,3) =0(1, 2) 5 (2, 3) g(l, 3) (23) 

Also, (17) is just the Ornstein-Zernicke equation in liq- 
uids 

K 2 f 

v{x) = v DH (x) + — — / dx v DH {x) v(x - x) (24) 
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(18) allows us to go beyond this equation systematically. 
The correlations introduced byK^O upset the superpo- 
sition principle (23) and the Ornstein-Zernicke equation 
(24). 



C. Virial Corrections to Linear Debye-Hiickel 
Screening 

At low density or small A we can systematically correct 
the linearized Debye-Hiickel limit. The purpose of it is to 
see how stable is the linearized theory for higher densities. 
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We note that an expansion in A of e AK in (18) is not 
an expansion in the coupling (although the diluteness 
is a property of the weakly coupled phase). Thus, we 
define the rcscaled and dimensionless variables using the 
short distance cutoff (core size) A = a 3 A, c = a 3 c and 
v = Vja? . The rescaled free energy is just 



F = 3£lnA - ^lnZ A 

v 

and the rescaled concentration per species is 
1 dlnZx 
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The first contribution in (25) is that of free particles while 
the second contribution is the excess/deficit free energy. 
The excess/deficit is caused by the interaction. 

The excess/deficit free energy can be virial expanded 
for low concentrations whatever the coupling, 
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in terms of which the concentration reads 
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The latter can be inverted to give 
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In terms of the concentrations, the free energy is then 
F = 3clnc - 3c + D 3 / 2 5 3/2 + D 2 c 2 + D b/2 c b/2 + ... (30) 
with 
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The virial expansion for arbitrary N c is involved. For 
N c = 2 all color integrations can be done explicitly. The 
result is 



and 



D 2 



£>3/2 = -^(E e «) 

24e ti Jo t 



(32) 



Q,/3 

e - 3£ (2-2e-3e 2 -9e 3 ) 
+e 3e (-2-2e + 3e 2 -9e 3 )) 



(33) 



and 



B 



54tt 3 / 2 



5/2 



\ a J \ a J 



with the classical binding energy 
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and e = ^/e a ep. In (34) only the logarithmic contribution 
was kept with the dimensionless cutoff A a reflecting on 
the infrared sensitivity of the virial coefficient. 

Z?3/2 is the linear Debye-Hiickel contribution to the ex- 
cess free energy while D 2 is the correction due to 2-body 
effect in the screened phase. The contributions e ±3e are 
the Boltzmann contributions from would-be molecules 
made of 2 charges a distant a apart. The virial expansion 
is exponentially sensitive to molecule formation, pointing 
to the quantum character of the sQGP even at low con- 
centrations, i.e. 
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For a 2-color sQGP, the Casimirs C 2q = C 2 q = 1/2 
and C 2g = 3 are of order 1 . Setting C 2 a s m 1 yields 
e i=a 1/aT. For a core a w 1/2T, the virial expansion 
around the linearized Debye-Hiickel limit breaks down 
for concentrations as low as c* w 10 _4 T 3 . 



IV. NONLINEAR POISSON-BOLTZMANN 
SCREENING 

To go beyond the linear Debye-Hiickel theory we need 
to solve the full mean-field equation (12), known as a 
Poisson-Boltzmann equation. In this section, we do so 
using the simplifying assumption that the non-linear so- 
lution remains spherically symmetric. This allows us to 
probe the stability of the linearized theory at strong cou- 
pling or higher densities. However, the reader should 
be prepared to see that at sufficiently strong coupling 
it would not be adequate either, with correlated charges 
forming non-spherical crystalline order. 

The non-linear radially symmetric Debye-Hiicfce^ equa- 
tion for the Abelian potential is 



4>"{r) +2<j)'{r)/r = 47r5n(0) 



(37) 



where the density variation in the r.h.s. can be written 
as a Boltzmann exponent* of the potential 



*The fact that we use the same F in both the Coulomb and 
core terms means that an increase in the coupling actually 
means a reduction of the temperature. 
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5n(4>) = n Lnt-Vcore) _ ^ ; (38 ) 

where the first term in the exponent is the Coulomb po- 
tential, defined to be positive. The second term, the core, 
is now needed to ensure that the density variation associ- 
ated with one particle is normalizable at small r. Below 
we set the normalization as 

J d 3 rSn((t)) = 1 (39) 
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FIG. 1. (Color online) (a) The function f(r) for 
r = 0.1,0.3,1,2,2.2 (top to bottom at the l.h.s.). (b) The 
integrand 4n r 2 5n(r) for the same values of the coupling. 
Larger couplings produce higher peaks. 



We set n = 1 below and seek a radial solution of the 
form 

0(0 = f(r) (40) 



with k = v47rT the Debye-Huckel inverse radius. One 
may think of f(r) as of an effective coupling. The 
value f(r) = 1 corresponds to the normalized solution 
of the linearized theory. If one includes the repulsive 
core at small distances, in weak coupling f(r) remains an 
r- independent constant, although different from 1. As we 
will see, at stronger coupling the deviations of f(r) from 
a constant value will display deviations from the linear 
Debye-Huckel theory. 

Inserting (40) into (37) allows for numerical solutions 
for weak and strong coupling. At large distances, the 
screened potential is weak and the linear Debye-Hiickel 
theory is valid with f(r) asymptotic a constant. The 
latter is fixed by the normalization (39). The numerical 
results for f(r) for different T are shown in Fig. 1(a). 
As expected, for T < 1 the function f(r) asymptotes a 
constant, except at small distances where the linearized 
theory never works. For T > 2 one enters a domain where 
/(r) is changing dramatically. We found that for T > 2 
the modification is fast growing, indicating a beginning 
of a phenomenon known as "over-screening" At large 
r our solution is large and oscillating, which indicate not 
only large deviations from a Debye theory but actually 
a complete breaking of the mean field approach itself. 
Such erratic solutions are a precursor for highly corre- 
lated state, beyond spherical mean field, and eventually 
a complete crystallization at very large T . 

The generalization of these results to a non-Abelian 
colored plasma is straightforward. In the SU(2) case, 
the original point charge can be thought of as having 
a particular color direction, while those in a screening 
cloud have any charges. As a result, there is an angular 
variable z ~ cos(6>) = Q ■ Qo in the Boltzmann factor 
as in (12). The potential &su2 averaged over the color 
orientations is 

^ = - 2 Jjzze — (41) 

where A = r<E> is the Coulomb part. The Casimir is ab- 
sorbed in r. These factors are reminiscent of the e ±3e 
in the virial expansion and have the same origin. (41) is 
plotted in Fig 2 along with its asymptote. The linear non- 
Abelian mean-field solution follows from the Abelian one 
by multiplying the squared Debye mass by < z 2 >= 1/3. 
Other gauge groups can be treated in the the same way. 



s After the work was basically completed the authors learned 
that it is widely used in important chemical and biological 
applications, see review in ref. [15]. 
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Indeed, in the SU(3) case the averaging over angles in- 
clude 6 variables (the 8 generators minus the 2 Casimirs) 
in the pertinent invariant measure. 




FIG. 2. (Color online) <&su2{r) (upper curve) with its small 
A/3 and large e A (A — 1)/2A 2 expansion. 



V. NONLINEAR SCREENING THROUGH A 
COULOMB CRYSTAL 

As we pointed out before, at very large T or very small 
T the system freezes into a crystal. To understand the 
difference between screening in the Debye-Hiickel limit 
and the screening in the crystalline phase, we discuss the 
induced potentials and the behavior of their pertinent 
partition functions. 

In the Debye limit the induced potential at small and 
intermediate coupling is of the form 

Vb(r) = -7 (l - ^ r/R °) « ~ (1 - r/2R D + ...) (42) 

The induced potential is finite at small r and rises lin- 
early with r. The pertinent classical partition function in 
the Einstein approximation (ignoring coupling between 
oscillations of different charges) is 

Z » J d 3 r e~ v ^ T » T 3 R ° (43) 

In a crystal the induced potential by all charges except 
one is quite different. If one starts with a pure Coulomb 
field in a cubic crystal, then 

"max / 1 Tlj 1 

V ^= E V -1 -T3 ( 44 ) 
* — ' r — n\ r 

n i — — n 7na:E 1 1 11 



where we introduced a cutoff of the sum n max and also 
subtracted the n x = n y = n z = term. The alternating 
signs cause the sum to converge. To smoothen-out the 
oscillations between even-odd finite n max we use instead 

T>/ \ V(n max ) + V(n max + 1) 

V[x, y, z) = (45) 

which is within few percent of the limit already at 
nmax = 3. The convergence is related to the fact that 
new distant charges added between n max and n max + 1 
do not have a nonzero dipole or quadrupole contribution. 
In fact they start contributing to the 1=4 octupole. As 
a result, the field is extremely flat at small r (near the 
origin) with only quartic r 4 corrections to the constant 
V(0), see Fig. 3. This is just Madelung constant which 
is the energy of the ideal zero T crystal (in our case it 
is 1.75). The figure shows also high multiples in the an- 
gular distribution, with alternating positive and negative 
(unstable) directions. 




FIG. 3. (Color online) Contour plot of the Coulomb poten- 
tial in the z=0 plane versus x and y (in units of the nearest 
neighbor distance) as induced by all charges except one in a 
cubic crystal. 

We note that the potential used in our molecular sim- 
ulations in [7] has in addition a repulsive core needed 
to prevent collapse at short distances. The result is a 
quadratic instead of quartic potential at the origin as is 
seen in Fig. 4, with V(r) w 9r 2 near the origin. The per- 
tinent partition function in the Einstein approximation 
is then 

Z» |d 3 re-^=(^) 3/2 , (46) 



G 



which leads to the usual classical potential energy of os- 
cillations 3T/2. However, as we have shown above, in 
order to understand the MD data one has to include the 
oscillations of the color vectors as well. 
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FIG. 4. (Color online) Induced potential with Coulomb 
and repulsive core in the z~0 plane: (a) along the crystal 
axes (green) and the diagonal (red), along with the fit (yel- 
low), (b) Planar contour plot for a cubic crystal. 



VI. FROM DEBYE-HUCKEL SCREENING TO 
DEBYE CRYSTAL 

As the density is increased or equivalently the coupling 
is increased or the temperature is decreased, the potential 
effects dominate over the kinetic effects in the cQGP. T 
which is the ratio of the potential to kinetic energy is then 
large. For simplicity and to compare with the molecular 
dynamics results in I, we specialize in this section to one 
species, say gluons only. Let T be the standard cQGP 
plasma parameter defined as [7] 
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with the Wigner-Seitz radius satisfying c(47ra^ s /3) 
1. In terms of (47) the excess free energy F cx is 



F(r) = f(o) + F cx (r) 



(48) 



and measures the Coulomb interacting part of the free 
energy. In the cQGP the excess free energy F ox per par- 
ticle is the Debye-Hiickel contribution 
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which is finite, attractive and nonlinear in the concentra- 
tion c. 

The higher order corrections in the concentration 
and/or coupling were discussed above and will not be 
repeated. Instead, here we note that (49) is finite and in- 
volves an integration over small as well as large momenta. 
As the density is increased in the cQGP there is increas- 
ing local ordering due to the strong colored Coulomb at- 
tractions causing the dilute and screened gas to undergo 
changes to a liquid and eventually to a crystal. 

Following Debye description of solids, we will assume 
that the allowed momenta in (49) are those commensu- 
rate with the total number of degrees of freedom, 
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which fixes the Debye cutoff qjj = 7r(9dc) 1 / 3 . For 
every particle moving in 3d space there are attached 
d = N c (N c — \)/2 independent but classical color charges. 
The cutoff Debye-Hiickel contribution is 
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The integrals are readily undone and the answer is 



(51) 
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with A — (2/3)(2/(7r d) 2 ) 1 / 3 . A similar construction was 
suggested by Brilliantov for the Abelian one component 
plasma [14]. The excess energy following from (52) is 
U ox = SFoH/cunr which is 



~(N 2 - 1) ^ ((ATf/ 2 atan(l/VAf )) (53) 
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For small T we recover the (linear) Debye-Hiickel limit, 
while for large T we obtain 
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where the Madelung constant is cmd ■ For N c 
A ps 0.392 and c MD = 2.64 



2: d 



Uex 




3/2 NT 

FIG. 5. (Color online) Potential energy as a function of 
kinetic energy, (3/2) NT for N=64 particles is N=64. The red 
points are the MD simulations, their statistical accuracy (not 
shown) is estimated to be 7 percents for all points. The blue 
line is a fit to eq. (55). the green curve is a fit to eq. (58). 



—0.65 + 2.25 T. 



(57) 



which is in in good agreement with the molecular dynam- 
ics result at small temperature or large T. In the dimen- 
sionless units of the Appendix, the conversion from T to 
T is given by T = 0.25/T. 

In the intermediate coupling regime with T ps 0.8 — 3, 
we enter the liquid regime of the cQGP (the right part 
of the figure as small T means high T) . Numerically, 



N 



-0.3- 



0.03 



T 



(58) 



The Debye-Huckel limit predicts V ex /N = -0.3/VT at 
large T or weak coupling but not the large negative con- 
stant. Its appearance can be traced back to the ap- 
pearance of a sharp peak in the inter-particle distance 
at about 1/n 1 / 3 as shown in Fig. lb in the liquid phase. 
In the gas phase, the the peak is spread in the form of a 
large Debye cloud. The negative constant contribution to 
the energy in (58) corresponds to the potential following 
from this peak. 



VII. SUMMARY 



This excess potential energy was calculated numeri- 
cally in I for the one-species cQGP with N c — 2 using 
molecular dynamics. Figure 5 shows the plot of the excess 
potential energy as a function of the total kinetic energy, 
(3/2)NT. The region of small temperature corresponds 
to large T while the region with large temperature to 
small r. For points corresponding to T ss 16 — 131 the 
excess energy can be approximated by ** 

re -0.5 + 2. 25T. (55) 

The numerical crystal analysis of section V yields the 
following excess energy per particle 

~ —0.44 + 2.5 T (56) 

with 3T/2 following from the oscillation energy in the 
harmonic potential, plus 1 T from the oscillations of the 
classical SU(2) color degrees of freedom ft. One can see 
that at least the specific heat of a low-T colored crystal 
is well reproduced. 

Using the definition of T in (47) in (54) the excess en- 
ergy per particle following from the Debye interpolation 
as a function of temperature is 



"See Appendix for the discussion of units. 
As explained in I, the unit color vector on a sphere is equiv- 
alent to an oscillator with one coordinate and one momentum. 



We have analyzed the screening and thermodynamics 
of a strongly coupled classical colored QGP. The basic 
question we asked in this paper was when and whether 
one can relate the textbook theories, such as Debye- 
Huckel screening or virial expansion, to our MD data. 

We have identified the onset of nonlinear screening us- 
ing the mean field approach, as well as onset of the cor- 
relations using a low density expansion around the mean 
field. We have carried explicitly the expression for the 
pressure of the cQGP to order c 5 / 2 . Not unexpectedly, 
we have found that the textbook methods are not really 
applicable in the liquid regime we studied. In particular, 
the spherically symmetric mean field treatment is very 
unstable already at medium coupling, and the virial ex- 
pansion fails at very small densities. Both are signals for 
the onset of a cluster formation (liquid) or long range 
ordering (crystal). 

We hope to address more issues in forthcoming publi- 
cations. In particular, issues related to the collective ex- 
citations of the system (sound and plasma-color waves), 
and propagation of external bodies in QGP (jet energy 
loss) etc. 

APPENDIX 

In this appendix we briefly discuss the units that were 
used in I: length, time and mass. The unit of length was 
chosen as the separation between two particles, r = A, 
at which the potential (the sum of Coulomb and core) 
is minimum. The density or concentration of particles 
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n = c, is dimcnsionless in these units. All simulations in 
I were carried with n = c = 0.3676. The unit of time To 
was set by the plasma frequency 

*>=f' = <£?>' /2 ' < 59 > 

with e 2 = Cia s in the non-Abelian case. The unit of 
mass was set by the particles energy m. 

All dimensional quantities can be expressed using these 
basic units. For instance, the kinetic energy is measured 
in units of toA 2 /tq . All simulations in I were carried with 
n = c = 0.3676 and A = 1, m = 1 and r = 1. With these 
conventions, the comparison with the molecular dynam- 
ics uses a ws = 0.866, C 2 a s = 0.216 and T = 0.25/T. 
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